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ABSTRACT 

We report the results of 19 years of Arecibo timing for two pulsars in the globular cluster NGC 5904 
(M5), PSR B1516+02A (M5A) and PSR B1516+02B (M5B). This has resulted in the measurement 
of the proper motions of these pulsars and, by extension, that of the cluster itself. M5B is a 7.95-ms 
pulsar in a binary system with a > 0.13 M companion and an orbital period of 6.86 days. In deep HST 
images, no optical counterpart is detected within ~ 2.5(7 of the position of the pulsar, implying that the 
companion is either a white dwarf or a low-mass main-sequence star. The eccentricity of the orbit (e = 
0.14) has allowed a measurement of the rate of advance of periastron: Co = (0.0142 ± 0.0007)°yr -1 . 
We argue that it is very likely that this periastron advance is due to the effects of general relativity, 
the total mass of the binary system then being (2.29 ± 0.17) M Q . The small measured mass function 
implies, in a statistical sense, that a very large fraction of this total mass is contained in the pulsar: 
M p = (2.08 ± 0.19) M e (1 a); there is a 5% probability that the mass of this object is < 1.72 M 
and a 0.77% probability that 1.2 Mq < M p < 1.44 Mq. Confirmation of the median mass for this 
neutron star would exclude most "soft" equations of state for dense neutron matter. Millisecond 
pulsars (MSPs) appear to have a much wider mass distribution than is found in double neutron star 
systems; about half of these objects are significantly more massive than 1.44 Mq. A possible cause is 
the much longer episode of mass accretion necessary to recycle a MSP, which in some cases corresponds 
to a much larger mass transfer. 

Subject headings: binaries: general — pulsars: general — pulsars: individual (PSR B1516+02A) — 
pulsars: individual (PSR B1516+02B) — neutron stars: general — equation of 
state: general 
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1. INTRODUCTION 

Over the past 21 years, more than 130 pulsars 
have been discovered in globular clusters (GCs) 5 . 
Among the first d iscoveries were PSR B 1516+02A and 
PSR B1516+02B (jWolszczan et al. 1989( 1. Both are lo- 
cated in the GC NGC 5904. This cluster is also known 
as M5 and we refer to these pulsars below as M5A and 
M5B. M5A is an isolated millisecond pulsar (MSP) with 
a spin period of 5.55 ms. M5B is a 7.95-ms pulsar in a 
binary system with a low-mass companion (see §5.2|) and 
an orbital period of 6.86 days. At the time of its dis- 
covery, this was the MSP with the most eccentric orbit 
known (e = 0.14), this being ~ 10 4 — 10 5 times larger 
than that of MSP- White Dwarf (WD) systems in the 
Galactic disk with similar orbital periods. 

In the Galactic disk, 80% of all known MSPs are found 
to be in binary systems and, with the single exception 
of PSR J1903+0327 (Champion et al. 2008), they are 
in low-eccentricity orbits with WD companions. In GCs, 
gravitational interactions with neighboring stars, or even 
exchange encou nters, can produce bin ary systems with 
eccentric orbits (jRasio fc Heggie 1995(1 . These high ec- 
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centricities allow the measurement of post-Keplerian pa- 
rameters such as the rate of advance of periastron (Co) 
or, in the future, the Einstein delay (7) that are not nor- 
mally measurable in MSP binaries in the Galactic disk. 
If these effects are relativistic, they allow estimates of the 
the total binary and component masses. 

When Anderson et al. (1997) published timing solu- 
tions for M5A and B they used the eccentricity of M5B 
to detect its periastron advance. However, the large rel- 
ative uncertainty of the measurement did not allow any 
astrophysically useful constraints on the total mass of 
the binary. In this paper we report the results of re- 
cent (2001 to 2008) 1.1-1.6 GHz (L-band) observations 
of M5. The first 2001 observations were part of an 
Arecib o search for pulsars in GCs, which found 11 new 
MSPs (|Hessels et al. 20071) . Three of these were found 
in M5, and subsequent observations of this GC were 
made chiefly with the aim of timing these new discov- 
eries (Stairs et al. 2008, in preparation). M5A and B are 
in the same radio beam as the new pulsars and they are 
clearly detectable in the L-band data, permitting timing 
(see Figs. Q]and[2]) of much better (M5A) or comparable 
(M5B) quality to that obtained at 430 MHz by Ander- 
son et al (1997). The whole dataset now spans nearly 19 
years and provides much improved timing parameters. 

2. OBSERVATIONS, DATA REDUCTION AND TIMING 

M5A and B were observed with the Arecibo 305-m 
radio telescope from the time of their discovery in 1989 
April until 1994 July using the 430-MHz Carriage House 
line feed. For this, a 10-MHz band centered at 430 MHz 
was used. The Arecibo correlation spectrometer made a 
3-level quantization of the signal and correlated this for 
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a total of 128 lags. These data were then integrated for 
506.58561 /is, and the orthogonal polarizations added in 
quadrature before being written to magnetic tape. The 
L-band observations began in 2001 June, using the "old" 
Gregorian L-Wide receiver {T sys = 40 K at 1400 MHz). 
The "new" L-Wide receiver (T sys = 25K at 1400 MHz) 
has been used since it was installed in the Gregorian 
dome in 2003 February. The Wide-band Arecibo Pulsar 
Processors (WAPPs, Dowd, Sisk & Hagen 2000) make a 
3-level digitization of the voltages over a 100-MHz band 
for both (linear) polarizations, autocorrelating these for 
a total of 256 lags. The data are then integrated for 
a total of 64 /is and the orthogonal polarizations added 
in quadrature and written to disk. At first, only one 
WAPP was available, and centered the observing band at 
1170 MHz or 1425 MHz. From 2003, three more WAPPs 
have been available, and we now use three of them to 
observe simultaneously at 1170, 1410 and 1510 MHz, the 
cleanest bands within the wide frequency coverage of the 
new L-Wide receiver. 

For all observations, the lags were Fourier transformed 
to generate power spectra. For the L-band observations, 
the power spectra were partially dedispersed at a disper- 
sion measure (DM) of 29.5 cm _3 pc and stored as a set of 
16 sub-bands on the disks of the Borg computer cluster 
at McGill University. At 1170 MHz, the partial dedis- 
persion introduces an extra smearing of 18 and 1.6/is for 
M5A and B respectively. Adding these values in quadra- 
ture to the dispersive smearing per channel (60.9 and 
59.7 /is respectively), we obtain a total dispersive smear- 
ing of 63.5 and 59.7 /is for M5A and B respectively, i.e., 
the sub-banding introduces very little extra smearing. 
The 430-MHz power spectra and L-band sub-bands were 
dedispersed at the known DM of these pulsars and folded 
modulo their spin periods. All the L-band data reported 
in this paper were processed using the PRESTO pulsar 
software package 6 . 

We added the best 1170 MHz detections of both pul- 
sars to derive "standard" pulse profiles, and these are 
displayed in Figure [U A minimal set of Gaussian curves 
were fitted to these profiles to derive synthetic templates 
for each pulsar, and these were then cross-correlated 
with each observat ion's pulse profile in the Fourier do- 
main (|Tavlor 19 92) to obtain topocentric times of arrival 
(TOAs). Adding all the 1170 MHz observations irre- 
spective of their signal-to-noise ratio (SNR, this varies 
from day to day because of diffractive interstellar scin- 
tillation), we obtain a "global" pulse profile which has 
a lower SNR than the "standard" pulse profile. We cal- 
culated the average flux densities for both pulsars (see 
Tabled]) from the off-pulse r.m.s. in their global profile, 
assuming a system equivalent flux density of 3.5 Jy which 
is valid for 1170 MHz at the high zenith angles required 
to observe M5. 

The TOAs were analyzed with TEMPO 7 , using the 
DE 405 Solar System ephemeris (|Standish 19981 ) to 
model the motion of the Arecibo radio telescope rel- 
ative to the Solar System Barycenter. The orbital 
parameters for M5B were m odeled using the Damour 
fe Deruelle orbital mod el (| Damour fe Deruelle 19851 
iDamour fe Deruelle 1986h . For most of the early 430- 

6 http:/ /www. cv.nrao.edu/~sransom/presto 

7 http: //www. atnf . csiro . au/research/pulsar/tempo/ 




Fig. 1. — Full-cycle pulse profiles for M5A (top) and M5B (bot- 
tom) at 1170 MHz, obtained by averaging the best detections of 
these pulsars at this frequency. The horizontal error bars denote 
the total time resolution of the system at 1170 MHz: 90.2 and 
87.5 fj.s for M5A and B respectively. 

MHz TOAs we have no reliable uncertainty estimates. 
Therefore, in order to achieve a reduced \ 2 of 1 for 
both pulsars, we attributed a constant uncertainty to the 
430-MHz TOAs that is similar to their unweighted rms. 
These times are listed in Table [T] With these TOAs, 
we obtain timing parameters that are virtually identical 
to those obtained by the previous analysis (Anderson et 
al. 1997). For the new L-band data, we find the TOA 
uncertainties to be under-estimated in the case of M5A. 
Multiplying these by a factor of about 1.5, we achieve 
a reduced x 2 of 1 for its L-band TOAs. In the case of 
M5B, this factor is 1.05. These values are similar to 
those derived for other MSPs timed with the same soft- 
ware (jFreire et al. 2008( 1. 

The resulting timing parameters and their l-er uncer- 
tainties are presented in Table [TJ We estimate these un- 
certainties to be twice the 1-a Monte-Carl o bootstrap 
(|Efron fe Tibshirani 1991 iPress et al. 1992f) uncertain- 
ties. We discuss the validity of this choice for the par- 
ticular case of the periastron advance of M5B in [J5] All 
the parameters that vary in time (a, 5, v and u>) are 
estimated for the arbitrary epoch, MJD = 54000 (2006 
September 22); To was the first periastron passage to 
occur after that date. The post-fit timing residuals are 
essentially featureless at the present timing accuracy (see 
Fig. [2]). The large gap without measurements between 
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TABLE 1 

Parameters for two pulsars in NGC 5904 



Observation and flux parameters 


PSR B1516- 


-02A 


PSR B1516+02B 


Start of 430-MHz observations 




47635 




End of 430-MHz observations 




49432 




Number of TOAs Q 430 MHz 


86 




82 


Residual rms & 430 MHz (/its) 


49 




114 


Start of L-band observations 




52087 




End of L-band observations 




54497 




Number of TOAs @ L-band 


1278 




162 


Uncertainty scale factor 3, 


1.50 




1.05 


Residual rms @ L-band (/is) 


9 




72 


Average flux density @ 1170 MHz (mjy) 


0.155 




0.027 



Ephemeris 



Reference Epoch (MJD) 


54000 


54000 


Right Ascension, a (J2000) 


15 h 18 m 33?32307(6) b 


15 h 18 m 31 ! r 4625(8) 


Declination, 8 (J2000) 


+02°05'27^' 435(3) 


+02° 05' 15'.' 30(3) 


Proper motion in a, fi a (mas yr _1 , J2000) .... 


4.6(4) 


3.4(1.2) 


Proper motion in <5, [ig (mas yr _1 , J2000) .... 


-8.9(1.0) 


-11.8(2.8) 


Spin frequency, u (Hz) 


180.063624055103(3) 


125.83458757935(6) 


Time derivative of u, v (10 — 15 Hz s — 1 ) 


-1.33874(4) 


0.05233(15) 


Dispersion Measure, DM (pccm -3 ) 


30.0545(10) 


29.46(3) 


Orbital period, Pi, (days) 




6.8584538(3) 


Projected size or orbit, x (1-s) 




3.04857(2) 


Orbital eccentricity, e 




0.137845(10) 


Time of passage through periastron, To (MJD) 




54004.02042(15) 


Longitude of periastron, oj(°) 




359.898(8) 


Rate of advance of periastron, Cj (° yr" 1 ) 




0.0142(7) 


Second time derivative of v, v (10 — 27 Hz s -2 ) 


[-0.6±0.5] c 


[4.4 ± 5.4] 


Time derivative of x, x (10~ 12 1-s/s) 




[-0.12 ± 0.09] 


Time derivative of P B , Pb (10~ 12 ) 




[15 ± 31] 


Derived parameters 


Spin period, P (ms) 


5.55359254401089(11) 


7.946940656275(4) 


Time derivative of P, P (10 -21 s s" 1 ) 


41.2899(13) 


-3.306(10) 


Mass function, /(M ) 




0.000646723(13) 


Total system mass, M (Mq) 




2.29(17) 


Maximum pulsar mass, M p max (Mq) 




2.52 


Minimum companion mass, M c m i n (Mq) 




0.13 



a This is the scale factor for TOA uncertainties required to achieve a reduced x 2 of unity in the 
solution. k 1-er uncertainties are presented in parenthesis. These are twice the estimate obtained using 
the Monte-Carlo bootstrap method. c Values in square brackets are not considered to be significant. 
They were not fit when determining the remaining timing parameters. 



Year 

1990 1995 2000 2005 




48000 50000 52000 54000 



Date (MJD) 

Fig. 2. — Timing residuals for M5A (upper) and M5B (lower), obtained with the ephemerides in Table[T] The large gap in time coverage 
is mostly due to the Arecibo upgrade. 
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the early 430-MHz data and later L-band data is in part 
due to the Arecibo upgrade of the late 1990's. We in- 
cluded an arbitrary time step between these two datasets 
in the fit. We have tested the timing solution by intro- 
ducing extra pulsar rotations between the two datasets, 
but these are always absorbed by this arbitrary time step, 
with no other changes in the fitted timing parameters. 

The positions, periods and period derivatives that we 
have obtained are consistent with those of Anderson et 
al (1997). In what follows, we analyze solely the newly 
measured parameters: the proper motions and the rate 
of advance of periastron of M5B. 

3. PROPER MOTIONS 

In the reference frame of a GC, the rms of the velocities 
of its pulsars along the orthogonal axes perpendicular to 
the line of sight should be the same as the rms of their 
velocities along the line of sight. The stellar rms velocity 
along the line o f sight at the center of M5 is 7.15 km s _1 
(iWebbink 19851) : the rms of the pulsar velocities should 
be smaller given the larger masses of the neutron star s 
(NSs). At the distance of M5, 7.5 kpc ([Harris 1996T ). 
the stellar rms velocity represents a proper motion rms 
of only 0.2 mas yr -1 . Given the present measurement 
precision, the estimated pulsar proper motions should be 
mutually consistent and reflect only the proper motion 
of the GC. The proper motion measurements in Table [1] 
are indeed ~ 1 a consistent with each other; from their 
weighted average we derive a proper motion for M5 of 
[i a — (4.3 ± 0.4) mas yr -1 and /is = (—9.6 ± 1.0) mas 
yr" 1 . 

M5 is one of the four GCs in the Galaxy for which opti- 
cal proper motion measurements have not provided con- 
sistent (i.e., agreeing within the formal uncertainty esti- 
mates) results (Dinescu, Girard & van Altena 1999). Our 
M5 proper motion measurement is in marginal agreement 
with the values derived by Scholz et al. [1996, \i a = 
(6.7 ± 0.5) mas yr" 1 and fi s = (-7.8 ± 0.4) mas yr -1 ], 
but is in good agreement with the values derived from 
Hipparcos [Odenkirchen et al. 1997, [i a — (3.3 ± 1.0) 
mas yr -1 and fig = (—10.1 ± 1.0) mas yr -1 ]. 

4. SEARCH FOR THE OPTICAL COUNTERPART OF M5B 

We have used the astrometric information on M5B 
to search for an optical counterpart in archival HST 
ACS/WFC data from programs GO 10120 and GO 
10615. The GO-10120 images were taken on 2004 August 
1 through the F435W, F625W and F658N filters. Since 
the uncertainty in the absolute astrometry of HST data 
is 1-2", we first tie the astrometry of the GO-10120 ACS 
images to the ICRS frame using UCAC2 stars (positional 
accuracy <0.070 " down to the magni tude limit of the 
UCAC2 catalog, IZacharias et al. 20041 ) . Since UCAC2 
standards in the small ACS field (3.4' x 3.4') are scarce, 
we use ground-based imaging of M5 to derive secondary 
standards. We retrieved from the public archive of the 
2.5 m Isaac Newton Telecope (INT) a 30-s Sloan-r Wide 
Field Camera image taken on 2004 June 8 and processed 
only the chip that contained the core of M5 (field of 
view ~23'xll'). Astrometric calibration of this image 
was achieved using 308 UCAC2 stars with positions cor- 
rected for proper motion to the epoch of the INT image. 
After fitting for shift, rotation angle, scale factor and dis- 
tortions, the final solution has rms residuals of 0.050" in 




Fig. 3.— A small portion of a GO-10120 HST image of M5, 
centered on the position of M5B, in negative. The picture was 
taken through the F435W filter of the Advanced Camera for Sur- 
veys/Wide Field Camera. The dark circle indicates the position of 
M5B and the radius (0.2") corresponds to the 2<r uncertainty in the 
absolute astrometry of the image. The nearest stars (white circles) 
lie at 2.7 and 3. Oct from the radio position; these are indicated as 
A and B and are discussed in J4] 

right ascension and 0.047" in declination. We selected a 
set of 198 secondary standards from unsaturated and rel- 
atively isolated stars in the INT image. These were used 
in turn to compute an astrometric solution for the short 
(70s) F435W distortion-corrected (using the multidrizzle 
software) exposure. The resulting fit for shift, rotation 
angle and scale factor has rms residuals of 0.017" in right 
ascension and 0.014" in declination. We estimate the fi- 
nal lcr accuracy of our ACS absolute astrometry as the 
quadratic sum of the errors in the UCAC2 astrometry, 
the UCAC2-INT tie and the INT-HST tie, i.e. 0.1" (or 
2 ACS pixels). 

The position of M5B at the epoch of the GO-10120 
observations is shown in Fig. [3l No optical sources are 
detected within the 2-a error circle of the radio po- 
sition (this only includes the uncertainty in the abso- 
lute astrometry as the uncertainties in the radio posi- 
tion are negligible). The photometry shows the near- 
est sources — indicated as "1" and "2" at distances 2.7 
and 3.0 er, respectively — to be main-sequence (MS) stars 
located -1.3 and ~6.1 mag (in F435W) below the MS 
turnoff. Using the M5 turnoff magnitude from Sandquist 
et al. (1996) and assuming B w FASbW, the 14-Gyr 
isochrones from Bergbusch & Vandenberg (1992) im- 
ply that "1" and "2" are ~0.75 and ~O.4M stars, re- 
spectively (for a distance modulus (m — M) — 14.41 
and E(B - V) = 0.03, see Sandquist et al 1996). To 
check for fainter stars, we stacked F435W images taken 
with fraction-pixel offsets from HST program GO 10615 
(2006 Feb 15) into a deep (8500 s) high-resolution (twice- 
oversampled) masterframe. Psf photometry revealed no 
additional stars within 3cr of M5B. A conservative upper 
limit for the detection limit is mF<i35W ~ 26 — 26.5 mag 
which corresponds to ^0.25-0.3M Q for MS stars in M5 
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(although close to star "1" the sensitivity is lower). 

In summary, we cannot exclude stars "1" or "2" as 
counterparts of M5B. However, the astrometry suggests 
that it is more likely that the true optical counterpart is 
fainter than our detection limit. That implies that the 
companion is either a WD or a faint, low-mass MS star. 

5. PERIASTRON ADVANCE OF M5B 

We have determined a highly significant estimate 
for the rate of advance of periastron of M5B: ill = 
0.0142(7)°yr~ 1 . The 1-er estimate provided directly by 
TEMPO is Co = 0.01422(43)°yr" 1 . To verify that these 
values are realistic, we have kept Co fixed, and fitted 
all the remaining timing parameters, recording the re- 
sulting x 2 - Doing this for a range of values of Co, we 
obtain the 1-er uncertainty as the half-width of the re- 
gion where [x 2 (w) — Y 2 (^min) < 11, Com\„ being the value 
that minimizes \ 2 (jSplaver et al. 20021 ). The result is 
Co = 0.001422(43)°yr~ 1 . Estimating all parameters us- 
ing a Monte-Carlo Bootstrap algorithm (Sj2]) we obtain 
Co = 0.01422(36)°yr- 1 . Furthermore, the lack of signifi- 
cant higher derivatives of the spin frequency (see Table[T]) 
suggests that, with the present timing precision, we are 
unable to detect any of this pulsar's timing noise which 
can contaminate the physical interpretation of its tim- 
ing uncertainties. This is a necessary pre-condition for 
an accurate estimation of parameter uncertainties. We 
therefore believe that the ler TEMPO uncertainty estimates 
are essentially accurate. We choose, however, to be more 
conservative by making our 1-cr uncertainties twice as 
large as the values suggested by the Monte-Carlo method 
(see Table [1} . This caution is due to our use of two dif- 
ferent datasets with a large time gap between them. 

This time gap is common to many Arecibo timing data 
sets, such as that of PSR J0751+1807. Nice et al. (2005) 
claimed a mass of 2.1 Mq for that pulsar. This claim has 
recently been retracted, the latest estimate for the pul- 
sar mass now being 1.261^2^0 ( lcr > Nice 2007 8 ). The 
problem with the earlier estimate was not related to the 
gap as the new estimate is based on essentially the same 
data. Its cause was the use of a necessarily imperfect 
ephemeris when folding the very earliest data. This re- 
sulted in orbital-dependent smearing of the pulse profiles 
that led to an error in the calculation of the orbital phase 
for the earliest data and an over-estimate of the orbital 
period decay. Because those data were folded online, 
this problem could only be solved by ignoring the earli- 
est TOAs. Our 430-MHz data contained the signals of 
more than one pulsar, so we had to record them to tape. 
This allowed us to re-fold them iteratively after the tim- 
ing solution had been obtained (Anderson et al. 1997). 
After updating the timing solution to 2008, we have also 
re-folded all our L-band data; none of the pulse profiles 
used in this work are therefore smeared due to imprecise 
folding. 

5.1. Is Co relativistic? 

The possible contributions to Co in a system contain- 
ing a pulsar and an extended star have been studied in 
detail by Lai, Bildsten & Kaspi (1995). In their analy- 
sis of the binary pulsar PSR J0045— 7319 they concluded 
that the only likely contribution to Co in such systems 

8 See also http://www.ns2007.org/talks/nice.pdf 



is from rotational deformation of the companion. If we 
assume that PSR J0045-7319 has a mass of 1.4 M , the 
companion mass is 8.8 Mq and its radius R c ~ 6.4 i? Q 
(jBell et al. 19951 ); the orbital separation a is ~ 126 R@. 
Using these values they reached the conclusion that the 
contribution to Co from tidal deformations, which are pro- 
portional to (R c /a) 3 ~ 1.32 x 10 -4 , is not significant in 
that system. For M5B, again assuming a pulsar mass of 
1.4 Mq with a MS companion of maximum mass 0.3 Mq 
(see §H]) and a radius R c = 0.3ro.3 Rq (where r , 3 ~ 1), 
then the orbital inclination i is ~ 24°, a ~ 18 Rq and 
{R c /a) 3 ~ 4.5 x 10 _6 r(}. 3 . This is 30 times smaller 
than for PSR J0045— 7319, becoming even less significant 
for smaller MS companion masses and correspondingly 
higher inclinations. 

Using equations 68 and 79 of Wex (1998) we can esti- 
mate the contribution to Co due to rotational deformation: 

kR 2 Cl 2 ( 3 

w ro t = n —T7 — - — ^tt ( 1 shi 2 9 + cot i sin 8 cos 9 cos <l>o 

ar(l — e z y \ 2 

(1) 

where n is the orbital angular frequency (2n/Pb = 1.06 x 
10 _5 rads ), k is the gyration radius (for a homogeneous 
sphere this is 0.63, while for any centrally condensed ob- 
jects this will always be smaller: it is about 0.2 for a 
completely convective star), f2 is the rotation rate rela- 
tive to break-up, 9 is the angle between the rotational 
and orbital angular momenta (if the companion is non- 
degenerate, these tend to be aligned, so 9 = 0) and <E>o is 
the longitude of the ascending node in a reference frame 
defined by the total angular momentum vector (see Fig. 9 
of Wex 1998). 

For the situation discussed above (a 1.4- Mq pul- 
sar with the largest possible MS companion, R c — 
O.3rrj.3-R0), we have R c /a = 16.55 x 10 _3 r .3. We can 
also calculate the break-up angular velocity: O max = 
y/GM c /R 3 = l/0.3y/GM Q /(r . 3 R & ) 3 = 2.09 x 
10~ 3 rQ~ 3 3 ^ 2 ' ) rads -1 . If the companion's rotation is 
tidally locked to its orbit around the pulsar, then D, = 
n/Q max = 5.06 x 10 _3 r^ 2 '. Therefore, Co mt = 2.65 x 
10 _15 r§ 3 rads _1 = 4.79 x 10 _6 (r^ 3 )°yr _1 (or three times 
this if the companion were to be a homogeneous sphere) . 
This is ~ 3 x 10 3 times smaller than the value we dedeter- 
mined (SJSJ). If the companion were to be significantly dis- 
tended for its mass (i.e., if r .3 > 1), as is the case for th e 
companion of PSR J1740-5340 (|Ferraro et al. 20011 ). 
then Co m t could be significant. This scenario can be ex- 
cluded, since no optical counterpart is readily detectable 
within ~ 2.5er of the pulsar (see SJ3]). 

If the companion were to be a WD, then it could be 
more massive than 0.3 Mq and still evade optical detec- 
tion. Irrespective of its mass, the contribution to Co from 
the tidal deformation of a WD is negligible, but that 
is not necessarily the case for the contribution from ro- 
tational deformation. As an example, we re-calculate 
Co mt for a 0.3-M w W D. For WDs, we have k = 0.45 
(|Livio fc Pringle 19 98) , more than twice as large as for 
fully convective stars. For WDs rrj.3 is of the order of 
0.1, i.e., the (R c /a) 2 term in eq. Q] would be ~ 10 2 times 
smaller than discussed for ro.3 ~ 1. However, a WD 
companion is not likely to be tidally locked. If it were 
spinning fast, then Cl ~ 1. While there is no special a 
priori reason why this should be true, it is a possibility 
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that cannot be excluded. This would mean that fl 2 could 
be ~ 4 x 10 4 times larger than discussed above, with Co TO t 
similar to the observed &). This is particularly so for the 
larger-sized WDs (those with the lowest masses). Fol- 
lowing Splaver et al. (2002), we note first that if the 
companion is not tidally locked, the angular momenta 
of the orbit and companion spin will probably not be 
aligned (6^0). In this case, the spin of the companion 
will induce a precession of the orbital plane. This will 
cause a change in i, affecting the projected semi-major 
axis of the orbit which will vary with a rate x. Rewriting 
equation 81 of Wex (1998), we can relate Co mt to x: 
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k'rot 



| sin 2 I 



— tanz- 
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This equation has the advantage that it does not depend 
on the mass (or the nature) of the companion. Thus, 
our observed 2-er upper limit of \x/x\ < 9.6 x 10 _14 s _1 
implies |ci ro t| < (1.7 X 10 _4 )°yr _1 times a geometric 
factor. In 80% of cases this geometric factor will be 
smaller than 10 and the upper limit for u> to t is similar 
to the present measurement uncertainty for Co. 

To summarize, u) I0 t can only be significant if the com- 
panion is degenerate, rotating near breakup velocity, and 
with its rotational angular momentum nearly aligned 
with the orbital angular momentum, making x unde- 
tectable. Otherwise Co is relativistic. 

5.2. Binary, pulsar and companion masses 

When Co is solely due to the effects of general relativity, 
we can measure the total mass of a binary system: 

1 

where T = GM /c 3 = 4.925490947^s. For M5B, we 
obtain M = (2.29 ± 0.17) Mq. For the nominal value 
of Co and a median i of 60°, the mass of the compan- 
ion is 0.173 M Q and the mass of the pulsar is 2.11 Mq. 
This is well above all NS masses that have been precisely 
measured to date. 

We calculated a 2-D probability distribution function 
(pdf) for the mass of the pulsar and the mass of the com- 
panion, assuming that the pdf for u> is a Gaussian with 
the half-width equal to the l-a uncertainty listed in Ta- 
ble [1] and an a priori constant probability for cosi. The 
two-dimensional pdf is then projected in both dimen- 
sions, resulting in 1-D pdfs for the mass of the pulsar and 
the mass of the companion. These are displayed graphi- 
cally in Fig. 31 The pulsar definitely has a mass smaller 
than 2.52 Mq, and the companion has a mass larger than 
0.13 Mq, the median and l-a limits for the pulsar and 
companion mass are 2mto.llMQ and Q.YI2^^Mq 
respectively. There is a 99%, 95% and 90% probabil- 
ity that the pulsar is more massive than 1.38, 1.72 and 
1.82 M Q respectively. There is a 0.77% probability that 
i is low enough to make the NS mass fall within the 
range of the components of double neutron star (DNS) 
systems: from 1.20 measured for th e companion 
of PSR J1756-2251 ([Faulkner et al. 2005h to 1.44 M^ 
measured for PSR B1913+16 (jWeisberg fc Taylor 20031 ). 
For M5B, assuming its nominal value of Co, these mass 
limits would imply that 7.9° < i < 10.2°. 



6. STATISTICAL EVALUATION OF MASS 
MEASUREMENTS 

M5B has the second largest mass estimate among all 
known NSs after PSR J1748-2021B (NGC 6440B). Be- 
cause of indeterminate orbital inclinations, all mass es- 
timates based solely on Co are probabilistic statements: 
one more PK parameter is necessary to have an unam- 
biguous determination of i and M p . No such parameters 
have yet been measured for M5B, NGC 6440B or any 
other eccentric MSP binaries in GCs; this is in some 
cases due to their low timing precision (like M5B, which 
is faint and has a broad pulse profile), and in others to 
their small timing baselines (like NGC 6440B). Neverthe- 
less, unambiguous upper limits for the pulsar masses and 
lower limits for the companion masses can always be ob- 
tained from a measurement of a relativistic Co alone (see 
Tabled]). Furthermore, in systems where the mass func- 
tion is very small and the total binary mass is very large 
(as for M5B and NGC 6440B) there is a much greater 
probability that most of the mass of the binary belongs 
to the pulsar itself, as described in ^5.21 

6.1. Evidence for high average neutron star masses 

An interesting feature of the eccentric binary MSPs in 
GCs is that as the binary mass increases, the mass func- 
tion / does not increase (see Table [5J . The exception 
is NGC 1851A; a system thought to have r esulted from 
an exchange interaction (|Freire et al. 20041 ). If these bi- 
naries were to have M p < 1.44 Mq then the increase in 
total mass would be due to higher companion masses, re- 
sulting in a general trend to higher mass functions. This 
is generally not the case. 

If we assume that all these GC MSPs have "normal" 
masses (between 1.2 and 1.44 Mq), we can then calculate 
the orbital inclinations of these binaries from their total 
masses and mass functions. These are displayed graphi- 
cally in Fig. [5] Of the five massive GC systems, four seem 
to have small orbital inclinations (i.e., with cosi > 0.8), 
the exception being NGC 1851A. A priori, one would 
expect only one out of five systems to have such a small 
i. We have used a Kolmogorov-Smirnov test to compare 
the cos i values corresponding to the nominal Co values 
and M p = 1.44M Q (only possible to calculate for the five 
massive binaries) with a fake set of 100 randomly ori- 
ented binary systems (i.e., with a uniform distribution of 
cos 2). We obtain a 0.46% probability that the observed 
distribution is extracted from this set with random or- 
bital orientations. If we use instead M p = 1.2Mq, we 
can calculate cosi for all the eccentric GC binaries. The 
probability that the resulting distribution of cosi is se- 
lected from the set with random orientations is 7.7 x 10 -5 . 

Such low inclinations might be less unlikely were there 
a tendency for pulsars to emit in a plane perpendicular 
to their orbit. During accretion, orbital angular momen- 
tum is transferred to the NS, making its rotation axis 
nearly perpendicular to its orbital plane. If the angle 
between the magnetic and rotational axes of pulsars (a) 
is small, the magnetic axis will describe a narrow cone 
nearly perpendicular to the orbital plane. However, no 
such tendency for small a has been described in the lit- 
erature. If a pulsar has a small a, its beam will proba- 
bly illuminate a smaller fraction of the sky, particularly 
if it is narrow. This can only make the low-a objects 
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Fig. 4. — Constraints on the masses of M5B and its companion. The hatched region is excluded by knowledge of the mass function and 
by sini < 1. The diagonal dashed lines correspond to a total system mass that causes a general-relativistic oj equal to, or within 1-cr of, 
the measured value. The six solid curves indicate constant inclinations. We also display the probability density function for the mass of 
the pulsar (top) and the mass of the companion (right), and mark the respective medians with vertical (horizontal) lines. For comparison, 
the gray bars indicate the range of "normal" NS masses (see £|5.2|1 . 



less likely to be detected 9 . Furthermore, no such ten- 
dency towards low orbital inclinations is seen among the 
two lighter binary MSPs in GCs, nor among the sys- 
tems with es timated orbital in clinations: NGC 6752A 
has i > 70° (iBassa et al._2006j); PSR J1909-3744 ha s 
i = 86.6(1)° (pacobv et al . 20051: iHotan et al. 2006T) : 
PSR J0437-4715 has i ~ 43° (|Verbiest et al. 20081) 
and the massive binar y PSR J1903+0327 has i ~ 79° 
dChamnion et al. 20081 ) 10 . 

9 The pulse profiles also tend to be wider in these nearly aligned 
rotators, further hindering their detection. We see no correlation 
between the "apparent" inclinations in Fig.[5]and the pulse- widths, 
although the pulsar profile of M5B is quite broad (Fig. [TJ. 

10 In many cases, high orbital inclinations enable the mea- 
surement of the mass of a MSP through Shapiro delay, produc- 
ing a selection of such inclinations among the systems with mea- 
sured mass. This should not be too strong in the sample we se- 
lected, systems with the timing precision of PSR J1909— 3744 and 



If the low mass functions of the massive binaries are 
not a result of systematically low orbital inclinations, 
they can only be due to systematically small companion 
masses. This implies that in the majority of these sys- 
tems the pulsar masses are significantly larger than in 
the lighter binaries. 

6.2. Implications for the equation of state of dense 
matter 

Because there is no physically plausible reason to as- 
sume that most massive binaries have a small i we now 
assume that the probability density of cos? is constant. 
We use this to calculate the mass pdfs from Co as de- 
scribed in § I5.2l for all the pulsars in Fig. [5] The pdfs are 
displayed graphically in Fig. [6] If the probability that 

PSR J0437— 4715 can have their Shapiro delays measured almost 
irrespective of their orbital inclination 
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TABLE 2 
Millisecond Pulsar Masses 



Name PSR 



GC 



P (ms) P b (days) 



//M© M/M a M c /M e M P /M Q Method b Ref.' 



Selected MSP Mass Measurements 



J07514 


-1807 




3.47877 


0.26314 


0.00000 


0.0009674 






i.26i; 2 


Pb, s 


1 


J1911- 


-5958A 


NGC 6752 


3.26619 


0.83711 


<0.00001 


0.002688 


l-5»_ . 10 


0.18(2) 




Opt. 


2 


J1909- 


-3744 




2.94711 


1.53345 


0.00000 


0.003122 


1.67^2 


0.2038(22) 


1.47±2 


r, s 


3 


J0437- 


-4715 




5.75745 


5.74105 


0.00002 


0.001243 


2.01(20) 


0.254(14) 


1.76(20) 


r, s 


1 


J19034 


-0327 




2.14991 


95.1741 


0.43668 


0.139607 


2.88(9) 


1.07(2) 


1.81(9) 


tj, s 


5 








Binary systems 


with indeterminate or 


bital inclinations 








J0024- 


-7204H 


47 Tucanae 


3.21034 


2.35770 


0.07056 


0.001927 


1.61(4) 


> 0.164 


< 1.52 


w 


6 


J1824- 


-2452C 


M28 


4.15828 


8.07781 


0.84704 


0.006553 


1.616(7) 


> 0.260 


< 1.367 


LU 


7 


J1748- 


-24461 


Terzan 5 


9.57019 


1.328 


0.428 


0.003658 


2.17(2) 


> 0.24 


< 1.96 


LD 


8 


J1748- 


-2446J d 


Terzan 5 


80.3379 


1.102 


0.350 


0.013066 


2.20(4) 


> 0.38 


< 1.96 


U) 


8 


B1516- 


402B 


M5 


7.94694 


6.85845 


0.13784 


0.000647 


2.29(17) 


> 0.13 


< 2.52 


UJ 




J0514- 


-4002A C 


NGC 1851 


4.99058 


18.7852 


0.88798 


0.145495 


2.453(14) 


> 0.96 


< 1.52 


Ul 


9 


J1748- 


-2021B 


NGC 6440 


16.76013 


20.5500 


0.57016 


0.000227 


2.91(25) 


> 0.11 


< 3.3 


dl 


10 



a Binary systems are sorted according to the total estimated mass M. b Methods are: Pj, - relativistic orbital decay, r, s - Shapiro delay 
"shape" and "range" , "Opt" - optically derived mass ratio, plus mass estimate based on spectrum of comp anio n, CJ - precession o f 
per iastron. c References are 1: dNice 20071) (total an d co mpanion masses n ot pr ovided), 2: (Bassa ct al. 2006), 3: (Hotan et al. 20061 ). 
4: l lVerbiest et al. 20081) . 5: llChampion et al. 2008D . 6: jFreire et al. 20031 ). 7: Begin et al. 20081) . 8: URansom et al. 20051) . 9: (Freire. 
Ransom & Gupta 2007); 10: JFreire ct al. 2008) d This pulsar is not technically a MSP, its spin period is longer than those found in most 
DNS systems. However, given the similarity of its orbital parameters to those of Terzan 5 I, we assume that it had a similar formation 
history. 6 Because of its large companion mass and eccentricity, this system is thought to have formed in an exchange interaction 
JFreire et al. 20041) . 
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Total Binary Mass (M Q ) 

Fig. 5. — Cosine of the orbital inclination i as a function of 
total binary mass for the eccentric MSP binaries in GCs. For each 
binary, the upper curve assumes a pulsar mass of 1.2 Mq and the 
lower curve 1.44 Mq. These are different for each binary because of 
their different mass functions. The vertical lines indicate the total 
masses corresponding to u) and its ±1,2 and 3 — a uncertainties, 
based on the most recent timing. If we assume randomly oriented 
orbits, then for any given total mass value the vertical distance 
between the two curves gives us the probability of the pulsar mass 
falling within the 1.2 — 1.44 Mq range. With the exception of 
NGC 1851A, these probabilities are significantly smaller for the 
more massive systems. 



1.2 Mq < M p < 1.44 Mq is small (as is the case for 
Terzan 5 I and J, M5B and NGC 6440B), it is a direct 
indication that the required orbital inclination ranges are 
very narrow (as shown in Fig. [5]) and therefore unlikely 
under the present assumption. For M5B, this probability 




1 1.5 2 2.5 3 
Pulsar Mass (Solar Masses) 

Fig. 6. — Probability distribution functions (pdfs) for the eccen- 
tric MSPs binaries in GCs (see also Fig. [5}. The mass pdfs of the 
MSPs in the less massive binaries (those with M < 2Mq ) are rep- 
resented by the dashed curves. The distribution of masses is much 
broader than is found for the components of DNS systems. Four 
out of a total of seven systems seem to be significantly more mas- 
sive than the most massive NS in DNS systems, PSR B1913+16. 

is only 0.77%, while for NGC 6440B it is even smaller, 
only 0.10%. Multiplying such probabilities for all the 
massive MSPs in GCs, we obtain a composite probabil- 
ity of 5.3 x 10~ 9 that all have masses between 1.2 and 
1.44 Mq. It is therefore very likely that some of these 
NSs are significantly more massive. 

As discussed above in H5.21 for M5B there is a 95% 
probability that the pulsar is more massive than 1.72 Mq. 
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This would exclude a third of the equations of state con- 
sidered in Lattimer & Prakash (2007). However, in this 
respect NGC 6440B should be far more constraining; 
there are 99 and 95% probabilities that the mass of that 
pulsar is > 2.01 and 2.36 Mq respectively. We consider 
the M5B result to be more secure: the non-detection of 
its companion (Q almost guarantees that its u> is purely 
relativistic ( jOj) . 

There are two MSPs listed in Table [2] for which large 
masses have already been determined, PSR J0437— 4715 
and PSR J1903+0327. The latter in particular has the 
potential for a precise, unambiguous measurement of a 
large pulsar mass in the very near future. At the mo- 
ment, it has not been confirmed whether its cl> is rela- 
tivistic or not, although it is likely to be so. These results 
strengthen the case for the existence of massive NSs. 

7. FORMATION OF MASSIVE NEUTRON STARS 

The MSP mass estimates in Tableland the mass pdfs 
in Fig. especially those of M5B and NGC 6440B, sug- 
gest that the distribution of MSP masses could span a 
factor of 2, a situation that is completely different to 
that found for the components of DNSs. NGC 1851A 
and the MSPs in the "light" (M < 2M ) binaries have 
masses smaller than 1.5 Mq, i.e., they are not signifi- 
cantly more massive than mildly recycled NSs, despite 
having spin frequencies of hundreds of Hz. In particular, 
the case of M28C demonstrated that if all NSs start with 
Mp > 1.2M Q , then some MSPs can be recycled by ac- 
creting < 0.15 M© from their companions. At the other 
end of the distribution, NGC 6440B could be twice as 
massive. 

It could be that MSPs were born with this wide 
range o f masses. Hydrodynamical core collapse si mu- 
lations (|Timmes et al. 19961: iBelczvnski et al. 20081) in- 
dicate that stars below ~ 18 Mq form ~ 1.20 — 1.35 M© 
NSs (such as 47 Tuc H, M28C and NGC 1851A), while 
stars with masses between 18 — 20 M Q form 1.8M Q NSs 
(similar to PSR J0437-4715, PSR J1903+0327, Terzan 
5 I, J and M5B). Above 20 Mq, stars experience par- 
tial fall-back of material immediately after the super- 
nova that can significantly increase the mass of the stel- 
lar remnant, making it either a super-massive NS (like 
NGC 6440B) or a black hole. 

This possibility raises the question of why such mas- 
sive NSs, while representing about half of the MSPs in 
Table [2j have not been found among the 9 known DNS 
systems. Most of the secondary NSs in DNS systems 
have masses between 1.2 and 1.3 Mq, and recently van 
den Heuvel (2007) suggested that these were formed by 
electron capture (EC) supernovae. The nine primary NSs 
in DNS systems were still likely formed in normal (iron 
core collapse) supernovae. The predicted percentage of 
massive NSs is quite small, and with only 9 known DNS 
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we are unlikely to see any massive NSs as the primary 
(jBelczvnski et al. 2008t ). a situation that is very different 
from what is derived for the MSPs. 

If the extra mass of some MSPs were instead acquired 
during the long accretion episodes that recycled them, 
we can explain naturally why we only see massive NSs 
as MSPs (not just in GCs, but also in the Galaxy, e.g. 
Verbiest et al. 2008, Champion et al. 2008) but not in 
DNS systems. 

8. CONCLUSION 

We have measured the positions and proper motions 
of M5A and B. This has allowed a detailed search for 
the companion of M5B. However, no object was detected 
within 2.5cr of the position of M5B to a magnitude limit 
of 26-26.5, indicating that its companion is either a low- 
mass MS star or a WD. We have measured the rate of 
advance of periastron for this binary system, concluding 
that it is very likely due solely to the effects of general 
relativity. In this case, the total mass of the binary is 
2.29 ± 0.17M Q , similar to the total masses of Terzan 5 I 
and J. Like those pulsars and NGC 6440B, the relatively 
low mass function for M5B indicates that most of the sys- 
tem mass is likely to be in the pulsar, which we estimate 
to be 2.08 ± 0.19 M (1 a). There is a 95% probability 
that the mass of this pulsar is above 1.72 Mq. If con- 
firmed, this would exclude about a third of the equations 
of state that are now accepted as possible descriptions of 
the bulk properties of super-dense matter. 

Together with other recent results, the large mass de- 
rived for M5B suggests that MSPs have a very broad 
mass distribution; half of these objects seem to be sig- 
nificantly more massive than 1.44 Mq. It is likely that 
all NSs began with the a narrow mass range like that 
found in DNS. They then accretted different amounts of 
matter (in some cases as much as their starting mass) 
during their evolution to the MSP phase. 
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